Quantum zigzag transition in ion chains 
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A string of trapped ions at zero temperature exhibits a structural phase transition to a zigzag 
structure, tuned by reducing the transverse trap potential or the interparticle distance. The transi- 
tion is driven by transverse, short wavelength vibrational modes. We argue that this is a quantum 
phase transition, which can be experimentally realized and probed. Indeed, by means of a mapping 
to the Ising model in a transverse field, we estimate the quantum critical point in terms of the 
system parameters, and find a finite, measurable deviation from the critical point predicted by the 
classical theory. A measurement procedure is suggested which can probe the effects of quantum 
fluctuations at criticality. These results can be extended to describe the transverse instability of 
ultracold polar molecules in a one dimensional optical lattice. 
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Crystals of singly-charged ions, ion Coulomb crystals, 
constitute a fascinating physical realization of the phe- 
nomenon discussed by Wigner for electrons in metals [l[ . 
These crystals result from the combined Coulomb repul- 
sive force and the restoring force of an external poten- 
tial, usually a harmonic trap Q • The strongly-correlated 
regime, leading to long-range order, is accessed by means 
of laser cooling, allowing one to obtain samples at tem- 
peratures T ranging from few mK to hundreds of fiK [3| . 

Ion chains have been most discussed in the literature 
focussed on quantum technological applications 0, 
They are realized in potentials with an enhanced elliptical 
geometry, so that the charges align along the trap axis 0- 
8]. Strictly speaking, in one dimension there is no crys- 
talline order even at ultralow temperatures. However the 
Coulomb interaction gives rise to very-slowly decaying 
density-density correlations, such that quasi long-range 
order is generally warranted . At a given interparticle 
distance the stability of the ion chain is determined by the 
strength of the transverse potential: if the corresponding 
frequency v t is below a critical value v c , a planar, zigzag 



structure is observed 
was shown that the classical structural transition linear- 
zigzag is of second-order. A Landau model was derived, 
proving that the soft mode of the phase transition is the 
zigzag mode of the linear chain (the transverse mode with 
the shortest wavelength) which drives the instability and 
determines the new structure (l3| . The corresponding 
Ginzburg-Landau equation in the continuum limit is re- 



111. In Ref. 12, hJ it 



ported in [14j . The results of these studies are strictly 



valid in the classical regime, when quantum and thermal 
fluctuations can be neglected. 

A question naturally emerges whether and how quan- 
tum fluctuations modify the behavior at the linear-zigzag 
instability. Theoretical studies on ion systems so far fo- 
cused on small samples [15]. In this Letter it is argued 
that the linear-zigzag transition is a quantum phase tran- 
sition, which can be mapped to the one of an Ising chain 
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FIG. 1: Phase diagram for a linear-zigzag transition, accord- 
ing to the mapping to the quantum Ising transition in ID, 
where T is the sample temperature and the dimensionless pa- 
rameter e is tuned by the confining potential or the interpar- 
ticle distance. The classical instability of the linear chain [13] 
occurs at e = 0. The quantum critical point, at e c > 0, 
separates the linear from the zigzag phase at T = 0. For 
< e < £ c quantum fluctuations dominate, and the crystal 
is in the linear (disordered) phase. The dashed lines indicate 
the boundaries of the quantum critical region. 



in a transverse field, describing a quantum ferromagnetic 
transition at T = [16[ . This equivalence was previously 
proposed in a theoretical study of linear-zigzag structures 
for electrons in quantum wires 



17] . Here we outline a 



mapping which allows us to estimate the quantum criti- 
cal point in terms of the system parameters, and deter- 
mine the experimental conditions required to observe it. 
The phase diagram is sketched in Fig.[TJ We find a finite, 
measurable deviation from the critical point predicted by 
the classical theory [13[. 

The system we consider is composed of N ions of 
mass m and charge Q. The ions mutually repel via 
the Coulomb interaction and are confined in the trans- 
verse plane (y — z) by an anisotropic harmonic poten- 
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tial, while periodic boundary conditions are assumed 
in the x direction. The ions experience the potential 
Vtot — V t + Vi n t > which is the sum of the transverse har- 
monic trap, V t — muf Y^,jLi Uj 1% an d the Coulomb re- 
pulsion, V int = (Q 2 /2)E l#J l/^(x l -x j ) 2 + (y l -y J ) 2 , 
while the motion in the z direction is frozen out by a 
tight confining potential, so that the system is essen- 
tially two dimensional. The parameters of the potential 
are such that the particles are aligned along the x-axis 
with uniform interparticle distance a, and their equilib- 



rium positions are x 



(0) 



,(°) 



0, with j = 1, 



N. 



"3 V 3. 

This configuration can be realized in a ring trap under 
the condition that one ion is pinned [fjj], or in a anhar- 
monic trap [lij . It also approximates the ions distri- 
bution at the center of a long chain in a linear Paul 

trap 0,11 HUES. 

At ultralow temperatures the ions perform harmonic 
vibrations about the mechanically stable equilibrium po- 
sitions. The linear chain is stable against longitudi- 
nal displacements from the equilibrium positions xf , 
while a soft mode in the transverse direction drives a 
structural transition into a zigzag configuration. This 
mode has wavelength A = 2a, such that in the clas- 
sical harmonic crystal the corresponding displacements 
are y^° ft — (— l^yo with yo the amplitude of the oscilla- 
tions. We focus on the dynamics of the transverse phonon 
modes of the linear chain close to the instability, and ex- 
pand V\ nt to fourth order in \yj\ -C a. Our starting point 
is the Landau model derived in Ref. [l3| ■ Close to the 
instability the displacement along y can be rewritten as 
y.j = (—l) J a4>j, where <pj is a dimensionless field which 
varies slowly with the ion position. To leading order in 
a gradient expansion, we obtain an effective potential for 
the field d>i 



We next define a parameter e 



2 )/u 2 . Classi- 



cally, the chain is stable for e < [13j. For e > 0, the 
local potential Vq in Eq. @ has minima at <f>± = ±<fi , 
with (fio = (y c j a)sjmej g\ the ions therefore order in a 
zigzag structure at either one of the configurations with 
= {—l) 3 4>±a, thereby breaking the reflection sym- 
metry at the a;-axis. This calculation, however, neglects 
quantum fluctuations, which induce tunnelling between 
the minima and tend to destroy the zigzag order. 

Quantum fluctuations are included in the model of 
the structural transition by means of a 1+1 dimensional 
quantum field theory [2(j. To this end, we write the par- 
tition function as 



Z = / V(j)e 



-S[<t>]/H 



(3) 



with the Euclidean action 



S[4>}= [ drJ2 
■Jo j=1 



N r 2 

ma 



K 



2 (d T ^y + v (cp J ) + -(h-cb J+1 y 



(4) 



where j3 = h/ksT — > oo for T — > 0. This model describes 
the quantum dynamics of the zigzag phonon mode in 
terms of a real continuous scalar field, — oo < 4>j{T) < oo. 
We now express it in a form which allows its mapping 
onto an effective model for the discrete field <Jj(r) = 
Sgn[</>j(r)]. We first divide the imaginary time (r) axis 
into M discrete steps separated by an infinitesimal in- 
terval of size St = /3/M. The partition function is then 
cast in the form Z = Tr {f M }, where T is a transfer ma- 
trix describing the time evolution from r to r' = r + St, 
whose elements read 
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where Vq is a local potential of the form 
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and the second term in Eq. (fT]) describes nearest-neighbor 
interaction at strength K = (Q 2 /a)ln2. The other pa- 
rameters are v c — <J 'AQ 2 / 'ma 3 and g — 24C$Q 2 / a 5 , 
with C a = J2j>i 1 /( 2 J + 1 ) a - Note that v c is the fre- 
quency at which the linear chain is unstable and a tran- 
sition to the zigzag configuration occurs within the clas- 
sical theory llj, ll3| . It is remarkable that the Coulomb 
repulsion is mapped to a short-range potential close to 
the linear-zigzag instability. This property is clearly 
visible in the spectrum of the transverse excitations: 



{K/m)k 2 [13(, characteristic of a 



chain with nearest-neighbor coupling and sound velocity 
c = Kjm. 



Here, T Joc (0 it St) = (<^| e - J W ? %-) is the propaga- 
tor describing the quantum dynamics of a particle sub- 
ject to the local Hamiltonian Hj = p 2 /(2m) + Vo((f)j), 
with Vb(</>) the double-well potential, Eq. @, and (j)j,Pj 
canonical conjugate variables. We now write the continu- 
ous field as 4>j = (f)Q<Tj +S(f>j , where <t| is the Pauli matrix, 
which can take values ±1 corresponding to the minima of 
the double-well potential, and S<f>j are the fluctuations. 
We then integrate over S<f>j , obtaining an effective trans- 
fer matrix in terms of Ising fields. When the system 
is well inside the zigzag phase, the typical frequency of 
tunnel-splitting in the double-wells, Aw, is much smaller 
than wq, the frequency of the oscillations about the well 
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minima. Under this approximation, we obtain 

TefAivjhWj}] (6) 
K^St 



N 



J^T/ oc [<r,-,<rJ]exp 



21) 



where the 2 x 2-matrix T; oc encodes the single particle 
dynamics within the double- well potential. The propaga- 
tor T toc obeys the symmetry Ti oc [a, a'] = T loc [—a, — a'], 
and can therefore be written as 



fi oc = const x (Aa° + Ba x 



(7) 



where o~ a are Pauli matrices in the basis where {a = ±1} 
denote the eigenvalues of a z , and <r° is the 2x2 unit ma- 
trix. The coefficients A and B acquire a simple form in 
the zigzag phase, when the ion wave functions are rela- 
tively well localized near ±</>o and can be approximated 
by Gaussians of width I. In this limit, using the inequal- 
ity e ~ 2 *o /i2 < tanh{Aw<5r/2} < 1, we obtain 



f loc ^T exp{Auja x ST/2} 



(8) 



with To a constant prefactor. As a result, 
the partition function reduces to the form Z « 
Zq J T>aexp (— Si[o~]/H), in which Zq is a constant and 
Si the action of an Ising chain in a transverse field 
subject to the Hamiltonian 



N 

i=i 
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■haf) 



(9) 



Here, the parameters are the fictitious exchange coupling 

J = K<{>1 = CjUps, (10) 

with Up = Q 2 /a the strength of the interaction and Cj = 
C3 In 2/(6(75), while the transverse field h = HAuj/2 is 
proportional to the splitting energy. Its dependence on 
the physical parameters can be estimated by means of a 
variational calculation [2(i | , which gives 



C h {UpU 2 k 



1/3 



(11) 



where Uk — h 2 /ma 2 is the typical kinetic energy scale 
of the atoms in the chain and Ch = (9C5/2) 1 / 3 . Equa- 
tion (TTTj) was obtained assuming Uk/Up ^> e 3 , which 
holds close to the critical points, where quantum fluctu- 
ations dominate. 

The model described by Hamiltonian © is known to 
exhibit a quantum phase transition at h/ J = 1 and T — 
0, separating an ordered phase at J > h (in our case, the 
stable zigzag configuration) from a disordered phase at 
J < h (the linear chain). In both phases the spectrum 
of excitations is characterized by a gap with energy A = 
2| J — h\. The corresponding phase diagram, showing the 



critical behavior as a function of s (and hence v t ) and 
T is sketched in Fig. Q] The critical regime (kpT > A) 
is characterized by universal power-law T-dependence of 
correlation functions with the critical exponents of the 
quantum transition [16[. 

The mapping to the Ising model, here presented, is 
strictly valid as long as the lowest energy levels of the 
local, double-well potential are below the barrier. In 
this regime, the tunneling between the wells is small and 
hence necessarily h <§; J, which is satisfied deep in the 
zigzag phase. The symmetry of the ordered phase was 
used to identify the effective model and estimate its pa- 
rameters. This procedure is strictly valid only in a higher 
dimension (4 — e). Nevertheless, from conformal symme- 
try arguments [2l| it is reasonable to assume that a quan- 
tum Ising model is the appropriate field theory for the 
l+l-D 4 -model describing the present system, and qual- 
itatively describes its critical behavior at h ~ J. This hy- 
pothesis is further supported by numerical studies [22j j . 
Under this assumption we estimate the quantum critical 
point in the regime of parameters where classical zigzag 
order is suppressed by quantum tunneling, and observe 
that it belongs to a universality class 16J, |2l| that differs 
from the one of the classical Landau model. 

The critical value of e, and thus of the transverse fre- 
quency v tl is dictated by the condition h = J, yielding 

n2/3 



C c {Uk/Up)*'* , C c = C h /Cj ~ 10 



(12) 



Hence, for given Uk/Up < 1, a quantum critical point 
is expected at e c >• Uk/Up. 

We now discuss what are the experimental conditions 
which are required in order to observe the quantum crit- 
ical point. Chains with dozens of ions are usually real- 
ized in linear Paul traps [D-IH , where the ion distribution 
along the chain is inhomogeneous. The considerations 
here reported apply in the center of the chain, when the 
chain contains several dozens of ions. More regular dis- 
tributions are expected in linear traps where the axial 
potential is anharmonic [18| . Here, one can tune through 
the critical point by either controlling the transverse fre- 
quency v t (typically in the MHz regime), or the spac- 
ing between neighboring ions a (typically several /im) by 
means of the axial confinement. To be able to distinguish 
the quantum disordered phase from the ordered (zigzag) 
phase, the frequency difference 5v = v c — v t should ex- 
ceed the experimentally accessible resolution. In order to 
estimate the required resolution, we write e c ~ 25v /v c in 
terms of a = ao x 1 um and of the ion mass m = nAm p , 
with Trip the proton mass and ua the atomic number, and 
find e c ~ 10 _4 /(n y iflo) 2/ ' 3 j leading to an upper bound on 
the frequency resolution 

Sum lQ-\u c /2){n A a )- 2 ' z (13) 

ranging from KHz (for protons with ao ~ 1) to several 
Hz for, e.g., 24 Mg + ions. This bound must also be com- 
pared with the heating time scale Th in ion traps, such 
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that SvTh <C 1 should be satisfied, which leads to de- 
manding conditions for the existing trapping setups [f|. 
A larger value 6u, and hence less restrictive conditions on 
observing the quantum critical point, could be reached 
in presence of screening, when the chain is embedded 
in a crystal of another ion species It is also essen- 
tial to reduce T below the energy scale characteristic of 
the gap away from criticality, i.e., A ~ h, resulting in 
T[mK] < W 3 C h (U$ c U P ) 1 / 3 /k B , which corresponds to 



T[mK] < 0.25 (n 2 A a%) 



-1/3 



(14) 



implying an upper bound of order ~ O.lmK for protons 
to several /iK for Mg + ions. 

The quantum critical point could be more easily ob- 
served in a quasi one-dimensional array of ultracold polar 
molecules, interacting via dipole-dipole repulsion, when 
the dipoles are aligned by an external field in a direc- 
tion perpendicular to the plane where their motion takes 



place [23|, |24J . Numerical simulations show that quan- 



tum tunneling clearly modifies the behavior at the linear- 
zigzag transition predicted by the classical theory [25j |. 
Assuming the dipoles are pinned by an optical lattice 
at fixed intcrparticle position, then the potential for the 
phononic modes at the instability is given by Eq. ([1} , and 
the mapping to the Ising model in the transverse field 
can be performed, where now Up = Cdd/(4ira 3 ), with 
Cdd the dipole interaction strength [25j, while C a ~ 10. 
Using typical numbers for the dipolar strength [23j one 
finds 5v ~ v c , showing that quantum fluctuations are 
significant at the linear-zigzag transition. 

The phase transition can be experimentally measured 
by light scattering via the structure factor S(k) [25J. In 
the critical region S(k) = Sf(k) + So5(k — ko), where 
feo = ir I a is the wave number of the ta%?,&%, Sq is pro- 
portional to the squared order parameter, So oc <pQ, while 
Sf is proportional to the isothermal susceptibility x 01 
the corresponding Ising system [l6j |. Approaching the 
quantum critical point in the zigzag phase (e > e c ) and 
for T < A, one has % ~ |e c - e|" 7 and $ ~ |e c - e\ 2f} 
where (3 — 1/8 and 7 = 7/4 are the exponents of the clas- 
sical two dimensional Ising model. In the critical region 
T > A (shaded area in Fig. 1), x ~ J 1 " 7 / 4 and cj>l ~ T 1 / 4 
Note that this critical behavior is strictly valid for 
an infinite system and assumes a uniform density of ions. 
It is expected to be approximately valid in an ion chain, 
when the variation of the density in a large region is neg- 
ligible within a correlation length. 

To conclude, we argue that the linear-zigzag instabil- 
ity in two-dimensional systems of trapped ions or polar 
molecules can be mapped to the one dimensional Ising 
model in a transverse field. This result demonstrates 
once more the potentialities offered by these systems as 
quantum simulators 0, [H, 0] , and more generally for 



quantum technological applications. 
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